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j^ ^ Abstract. We present a global analysis of available data on inclusive structure 

functions measured in electron-proton scattering at small values of Bjorken-x, including 
the latest data from the combined HERA analysis on reduced cross sections. Further, 
we discuss the kinematical domain where significant deviations from NLO-DGLAP 

^Sl should be expected and the ability of non-linnear physics to account for such deviations. 

r^ In the limit of small Bjorken-x, deviations from standard coUinear perturbation 

O theory are expected on account of large gluon densities. Therefore, non-linearities 

,-H become relevant and need to be accounted for in the evolution. A solid theoretical 

>- framework in which to address the non-linear dynamics driving the small-x evolution of 

K> the proton (nucleus) wave function is provided by the Colour Glass Condensate (CGC) 

H [H El |3], along with its B-JIMWLK evolution equations (HE]. On the phenomenological 

side, the implementation of saturation effects has proceeded via dipole-models in which 

the general physical properties of the CGC are implemented in an effective manner. 

The Balitsky-Kovchegov equation [5l [6], simplest theoretical realization of the 
CGC, fails to reproduce experimental data. However, it was demonstrated in ^ that 
considering only running coupling corrections accounts for most of the NLO effects. The 
ability of the rcBK equation to correctly describe experimental data was shown in [S] . 
The rcBK equation, the most reliable and phenomenologically usable non-linear small-x 
evolution tool, has become the standard choice for CGC computation of experimental 
observables (e.g., ^[101111]). 

In the framework above described, a global fit of the inclusive DIS structure function 
led to the public release of a parametrization of the dipole scattering cross section at 
small values of x [HI [12], which we now proceed to present. 



AAMQS setup 

The reduced DIS cross-section a^ can be written in terms of the virtual photon-proton 
cross-section aT,L, being y = Q'^/sx the inelasticity variable, as follows 

In the dipole formulation of QCD, valid for small-x, the photon-proton cross-sections 
for transverse and longitudinal polarization of the virtual photon can be expressed as 

r-l 

where "i/xi is the light-cone wave function for a virtual photon to fluctuate into a qq 
dipole of quark flavor / (with mass ruf and electric charge e/). A/'(r, x) is the imaginary 
part of the dipole-target scattering amplitude averaged over impact parameter, with r 
the transverse dipole size. The transverse area over which quarks of a given flavour 
are distributed is represented by (twice) ctq,/. Light quarks are taken to be identically 
distributed {(Joj=u,d,s = o'o,iight), while when accounting for heavy flavour contributions 
in ([I]) we allow aQj=c,b = <^Q,heavy to be different from ao,Hght- (see [12] for a more detailed 
explanation). The quantities ao^ught and ao^heavy are free fit parameters. 
The evolution of A/'(r, x) is given by the rcBK equation: 

dAf{r, x) 



|rfriir™-(r,ri,r2) 



d ln(xo/x) 

X [Af{ri, x) + M{r2, x) - Af{r, x) - Af{ri,x) X{r2,x)] , (3) 

with K^^"^ the evolution kernel including running coupling corrections [13j . The running 
coupling in K™^ in ^ is evaluated at 1-loop accuracy in coordinate space 

, 2s 47r ^ 2 , . 

yy^ ) = 7, ,__ / 4^72 \ ' /^0,n^ = 11 - -Uf , (4) 
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where the constant C accounts for the uncertainty in the Fourier transform from 
momentum to coordinate space and is also a fit parameters. The coupling is evaluated 
with Uf = 3 when only light quark contributions are taken into account, and with a 
variable flavour scheme once heavy flavours are included. Since in the rcBK equation 
(pj) all dipoles sizes are explored, an infrared regulation is called for: ^^(r^ > rj^) = 
afr = 0.7 where rj-^ is the dipole size at which the coupling reaches afr- 

We have explored in [12] different families of initial conditions to solve (IS]), finding 
that the results show a negligible dependence on the specific choice. For brevity, only 
results obtained with the GBW initial condition 



M ^(r, a: = Xo) = 1 — exp 



(5) 
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are shown here. The fit parameters QIq and 7 are, respectively, the saturation scale 
and the characteristic fall-off of the dipole scattering amplitude with decreasing r. 



Results 

In Fig. fl] (left) we show the comparison of data for reduced cross section with the 
results of a fit to the combined Hl/ZEUS data, E665 and NMC data for F2, and 
available data on the charm contribution to F2 and a^. within cuts x < 0.01 and 
Q^ < 50 GeV^. The parameters resulting from such fit were: x^/d-o.f. = 1.30, 
Qlo,H,htihea.y) = 0.23(0.22)GeV2, ao(.) = 36.36(20.38)mb, 7(0 = 1.24(0.92), C=7.86. 
Notice the remarkable agreement of the AAMQS calculation with the experimental data 
despite their high accuracy (errors are often smaller than the size of the symbols used 
for the data points). In Fig. [I] (right) we present a comparison to the charm component 
of F2^c and ar,c with the results of the fit described above, showing a good description 
of the experimental data as well. 





Fit including lieavy quarks 
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Comparison with data on F and o^. 
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Figure 1. (left) comparison of ct^ experimental data (black squares) with our fit 
results (red circles), (right) Coparison of experimental data for i^2c (black squares), 
and arc (red squares) with our resuts (cyan circles) 



As an independent check, since no data of this sort are included in any of the fits, 
we calculate the longitudinal structure function Fl and compare our result to the latest 
HI measurement [13]. The result is shown in Fig. ^as a function of Q^ for the different 
values of x finding an excellent agreement. 

Finally, we discuss the kinematical domain where significant deviations from NLO- 
DGLAP should be expected. In [15] it was shown that there is tension in NLO-DGLAP 
fits to HERA data on structure functions when data sets with Q^ below the estimated 
saturation scale of the proton were excluded from the analysis. Within the AAMQS 
framework a similar exercise can be performed by excluding from the fits all data above 
a given x = Xcut, and then extrapolating to the unfitted region. Since linear and non- 
linear effects are expected to be present in different (and complementary) regions of 
phase space, we need to determine where the saturation boundary is. To do so, in [16] 
both linear and non-linear fits to partial sets of data are performed in their respective 
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Figure 2. Comparison of the 
latest Fl data [Tl] (red dots) with 
the AAMQS calculation using the 
results of the fit shown if Fig. fll (blue 
triangles). 



safe (x, Q^) regions, leaving a common unfitted region in which the extrapolation from 
both formalisms is compared. We find that the deviations from NLO-DGLAP reported 
in [15] do not appear when the non-linear rcBK fit is performed at low-x and extrapolated 
to the unfitted region. This suggests the presence of saturation effects which are not 
accounted for in linear evolution. 
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